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ABSTRACT 

The seemingly unrelated regression (SUR) proposed by Zellner consists ofL regression equations each of which 
satisfies the assumptions of the standard regression model. These assumptions are not always satisfied mostly in 
Economics, Social Sciences and Agricultural Economics which may lead to adverse consequences on the 
estimator parameters properties. Literature has revealed that multicollinearity often affects the efficiency of 
SUR estimators and the efficiency in the SUR formulation increases, the more the correlations between error 
vector differ from zero and the closer the explanatory variables for each response being uncorrelated. This 
study therefore examined the effect of correlation between the error terms and autocorrelation on seven 
methods of parameter estimation in SUR model using Monte Carlo approach. 

A two equation model was considered in which the first equation has the presence of autocorrelation 
and correlation between the error terms exists between the two equations. The levels of correlation between the 
error terms were specified as CR = -0.99, -0.9, -0.8, -0.6, -0.4, -0.2, 0, +0.2, +0.4, +0.6, +0.8, +0.9 and +0.99 
and autocorrelation levels RE= -0.99, -0.9, -0.8, -0.6, -0.4, -0.2, 0, +0.2, +0.4, +0.6, +0.8, +0.9 and +0.99 A 
Monte Carlo experiment of 1000 trials was carried out at five sample sizes 20, 30, 50, 100 and 250. The seven 
estimation methods; Ordinary Least Squares (OLS), Cochran - Orcut (GLS2), Maximum Likelihood Estimator 
(MLE), Multivariate Regression, Full Information Maximum Likelihood (FIML), Seemingly Unrelated 
Regression (SUR) Model and Three Stage Least Squares (3SLS) were used and their performances were 
critically examined. Finite properties of estimators' criteria such as bias, absolute bias, variance and mean 
squared error were used for methods comparison. 

The results show that the performances of the estimators cannot be solely determined by the evaluation 
given by bias criterion because it always behaves differently from other criteria. For the eight different cases 
considered in this study, it was observed that when the sample size is small (i.e 20 or 30) there is high 
variability among the estimators but as the sample size increases the variance of the estimator decreases and 
the performances of the estimators become asymptotically the same. 

In the presence of correlation between the error terms and autocorrelation, the estimator of MLE is 
preferred to estimate all the parameters of the model at all the level of sample sizes. 
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I. INTRODUCTION 

The seemingly unrelated regression (SUR) model is well known in the Econometric literature (Zellner, 1962, 
Srivastava and Giles, (1987), Greene (1993) but is less known elsewhere, its benefits have been explored by several authors 
and more recently the SUR model is being applied in Agricultural Economics (O' Dorell et al 1999), Wilde et al (1999). Its 
application in the natural and medical sciences is likely to increase once scientists in the disciplines are exposed to its 
potential. 

The SUR estimation procedures which enable anefficient joint estimation of all the regression parameters was first 
reported byZelfner (1962) which involves the application of Aitken's GeneralisedLeastsquares(AGLS), (Powell 1965) to the 
whole system of equations.Zellner (1962 & 1963), Zellner&Theil (1962) submitted that the joint estimationprocedure of 
SUR is more efficient than the equation-by-equation estimationprocedure of the Ordinary Least Square (OLS) and the gain 
in efficiency would bemagnified if the contemporaneous correlation between each pair of thedisturbances in the SUR system 
of equations is very high and explanatoryvariables (covariates) in different equations are uncorrelated. In other words, 
theefficiency in the SUR formulation increases the more the correlation betweenerror vector differs from zero and the closer 
the explanatory variables for eachresponse are to being uncorrelated. 



THE IJES 



www.theijes.com 



The IJES 



Page 27 



Autocorrelation effects of seemingly unrelated regression (sur) model on some estimators 



After the much celebrated Zellner's joint generalized least squares estimator, several other estimators 
for different SUR systems were developed by manyscholars to address different situations being investigated. 
For instance, Jackson(2002) developed an estimator for SUR system that could be used to modelelection returns 
in a multiparty election. Sparks (2004) developed a SURprocedure that is applicable to environmental situations 
especially when missingand censored data are inevitable. In share equation systems with randomcoefficients, 
Mandy & Martins-Filho (1993) proposed a consistent andasymptotically efficient estimator for SUR systems 
that have additiveheteroscedastic contemporaneous correlation. They followed Amemiya (1977) byusing 
Generalized Least Squares (GLS) to estimate the parameters of the covariance matrix. Furthermore, Lang, 
Adebayo &Fahrmeir (2002), Adebayo(2003), and Lang et al (2003) in their works also extended the usual 
parametricSUR model to Semiparametric SUR (SSUR) and Geoadditive SUR models withina Bayesian context. 
Also O'Donnell et al (1999) and Wilde et al (1999) developedSUR estimators that are applicable in Agricultural 
Economics. More recently, Foschi (2004) provided some new numerical procedures that could successively and 
efficientlysolve a large scale of SUR model. .In all the estimation procedures developed for different SUR 
situations asreported above, Zellner's basic recommendation for high contemporaneouscorrelation between the 
error vectors with uncorrelated explanatory variableswithin each response equations was also maintained. 
However, in most practical situations, the explanatoryvariables across the different equations in SUR systems 
are often correlated. Also, it may be necessary to jointly regress the demand for two or morecomplementary 
products like automobiles and gasoline on peoples' income andexpenditures on other products within the SUR 
framework. While the twodemands (responses) would obviously correlate through their error, satisfying thefirst 
basic requirement of SUR estimation, people's income and their expenditureon other products should not be 
expected to be uncorrelated thereby, violatingthe second important condition. Therefore, the existence of this 
kind ofrelationship needed to be recognized and accorded proper management withinthe SUR context such that 
the efficiency of SUR estimator would not becompromised.lt is now obvious, due to several instances of SUR 
highlighted above, that theindependent variables are often correlated (collinear). 

The seemingly unrelated regression proposed by Zellner (Zellner; 1962) consists of L regression 
equations each of which satisfies the assumptions of the standard regression model: 

. . (1.1) 

y L =x l Pl+Ul 

Wherey ; and U; are Nxl vectors and X; is aNxK matrix. Notice that by stacking all equations together we can 
write this system as 



>'i 



0. 



.0 



0 



.X 



p, 



\HLJ 



LNxl LNxXiki likixl LNxl 

Y = Xp + U 



y u Lj 



(1.2) 



(1.3) 



So we stack all the equations together into a system of the OLS form in (1.3). This suggests as an estimation 
procedure to run OLS on this system, i.e to consider 
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1.2 Efficiency of SUR model for estimating regression coefficients. 

Note that if either (7 = 0 for all j ^ i or X ; = X ^ for all j ^ i , then, the two model formulations produce 
estimator identical to 



0j=( X J X j)' X 'j Y J 



forj=l,2,...,p 



The efficiency in the SUR formulation increases, the more the correlations between error vector differ 
from zero and the closer the explanatory variables for each response being uncorrelated (e.g. see Sparks (1987) ) 
discussed how to select variables and parameter estimators for the SUR model. The standard errors for the set of 
the regression parameter estimates in the SUR formulation are given by the diagonal elements of (X'^ 'X) 1 

while for the unrelated formulation they are the appropriate diagonal elements of (XjX) 1 for j=l,2...,p. 

These can be used to gain an idea of the relative merits of the SUR model formulation for estimating the 
regression parameters of the model. 

Generally when CT.. = 0 and the covariances are known, it can be shown that the diagonal elements 

of (X'jX) 1 are larger than the corresponding diagonal elements of (X'£ "'X) _1 for each j. The mean square error 

of prediction using the generalized least squares estimate is smaller. This is not generally true when the 
covariances are unknown but depend on the sample size n (Zellner, 1963; Kmenta and Gilbert, 1968; Revankar, 
1974 ; Mehta and Swamy ,1976; and Maeshiro ,1980). When fitting regression models with small sample sizes, 
is unlikely that the seemingly unrelated regression formulation and related generalized least squares estimates 
are going to add much value. 

However, the larger the sample size, the more reliable the estimate of cr and hence the more likely an 
advantage is gained from the seemingly unrelated regression formulation of the model. 



Consequently, this study examines the performances of some estimation methods in Seemingly Unrelated 
Regression model in the presence of autocorrelation with the intention of studying their effects on the estimators 
and identifying the preferred estimator(s) of the model parameters. 

Very specifically, the study aims at the following: 

(i) Examine the effect of sample size on the performance of the estimators 

(ii) Examine the effect of autocorrelation (RE) and correlation between the error terms (CR) on the 

performance of seven estimators. 

(iii) Identify the estimator that yields the most preferred estimates under separate or joint influence of the 

three correlation effects under consideration. 

II. THE MODEL FORMULATION 

The Seemingly Unrelated regression (SUR) Model used in this research work is given as 

y it = Ah + PuXu + Pn x it + u i< (3-1) 
where U h = pli X{t _ X) +e lt , e lt «(0,(T 2 ). 

y 2 , = Pai + Pn x u +A 2 *3, + u 2, > u 2, ~ N(0, cr 2 ) (3.2) 

NOTE: (1) Autocorrelation exists in equations (3.1) 

(2) There is correlation between Ui and U2 of the two equations 
(3) There is no correlation between Xj and X 3 in equation (3.2), thus, equation (3.2) appears 
as control equation. 
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III. EQUATION USED FOR GENERATING VALUES IN SIMULATION 

The equation used for generating values of the variables in the simulation study as proposed by Ayinde 
K.(2007) is given below 

Suppose W t ~ N(jU,<jf) i = 1,2 . If these variables are correlated, then, W, and W 2 can be generated with 
the equations 

I 2 (33) 

W 2 =/u 2 +pa 2 z l +a 2 z 2 ^l-p 2 

whereZ;~N(0,l) i = 1,2 and \p\ < 1 is the value of the correlation between the two variables. 

3.1 Other Specifications 

1. Sample Size(n) of 20, 30, 50, 100 and 250 were used in the simulation 

2. The following levels were used for the correlations studied: 

a. Autocorrelation(RE) : -0.99, -0.9,-0.8, -0.6, -0.4, -0.2, 0, 0.2, 0.4, 0.6, 0.8,0.9, 0.99 

b. Correlation between error term (CR) : -0.99, -0.9,-0.8, -0.6, -0.4, -0.2, 0, 0.2, 0.4, 0.6, 0.8, 0.9, 0.99 

c. Replication (RR) : we make use of 1000 replications 

d. Two RUNS were done for the simulations which were averaged at analysis stage. 

3.2 Criteria for comparism 

Evaluation and comparison of the seven (7) estimators were examined using the finite sampling properties of 
estimators which include Bias (BB), Absolute Bias (AB), Variance (VB) and the Mean Square Error (MB) 
criteria. 

A 

Mathematically, for any estimator J3 f of /? of model (3.1) & (3.2) 



A 1 7? A 

R M 



(ii)Bias fi t =-£ /?,-/?, =J3-J3 i 

\ J K y=l V J 



(iii)AB B t =-£ 

7=1 



v j R 



P,-Pj 



(iv)Var 



f A \ 
Pi 

V J 



1 ^ 



R 



7=1 



A A 

fit,- Pi 



(v) mse\ P^\ = ^Y\PirPi 



J R 7=1 



, fori = 0, 1, 2 andj = 1,2,... , R. 



Using a computer program which was written with TSP software package to estimate all the model 
parameters and the criteria, the performances of seven estimation methods; Ordinary Least Squares (OLS), 
Cochran - Orcut (COCR), Maximum Likelihood Estimator (MLE), Multivariate Regression, Full Information 
Maximum Likelihood (FIML), Seemingly Unrelated Regression (SUR) and Three Stage Least Squares (3SLS) 
were examined by subjecting the results obtained from each finite properties of the estimators into a multi factor 
analysis of variance model. Consequently, the highest order significant interaction effect which has "method" as 
a factor is further examined using the Least Significance Difference (LSD) test. The estimated marginal mean of 
the factor was investigated out at a particular combination of levels of the correlations in which estimators were 
preferred. An estimator is most preferred at a particular combination of levels of the correlation if the marginal 
means is the smallest. All estimators whose estimated marginal means are not significantly different from the 
most preferred are also preferred. 



IV. RESULTS 

The summary of results from the Analysis of variance tables of the criteria showing the effect of the 
estimators, correlation between the error term sand autocorrelation on P; are presented in Table 4.5.1 below. 
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TABLE 4.5.1: ANOVA for sample size of 20 



n 


sov 


EQN 


Pi 


df 


TYPE III SUM OF SQUARES 












Bias 


Absolute Bias 


Variance 


Mean Square 


20 


RE 


1 


(301 


12 


892.446*** 


115926.509*** 


2445822.237*** 


3951716.298*** 








(311 


12 


029*** 


32.515*** 


95.927*** 


96.084*** 








(321 


12 


.012 


24.373*** 


23.509*** 


87.017*** 






2 


(302 


12 


.112 


103.206*** 


122116.658*** 


128548.527*** 








P12 


12 


.063 


.628*** 


093*** 


091*** 








(322 


12 


232*** 


.605*** 


1 13*** 


125*** 




CR 


1 


(301 


12 


.670 


.003 


.003 


.005 








(311 


12 


.001 


6.016E-5 


8.532E-5 


8.897E-5 








P21 


12 


7.468*** 


3.807*** 


3.176*** 


5.004*** 






2 


(302 


12 


3.519 


45.130*** 


113879.706*** 


119769.347*** 








(312 


12 


.513*** 


.032 


.224*** 


032*** 








(322 


12 


3.006*** 


404*** 


139*** 


on *** 




M 


1 


(301 


6 


315.786*** 


83483.317*** 


4080093.223*** 


6466311.896*** 








(311 


6 


.000 


4.612*** 


5 977*** 


5.990*** 








P21 


6 


.007 


5.564*** 


2.320*** 


9.779*** 






2 


(302 


6 


.042 


45.091*** 


232859.705*** 


243905.100*** 








P12 


6 


.476*** 


141*** 


002*** 


.006*** 








P22 


6 


.086*** 


2.096*** 


.361*** 


391 *** 




RE*CR 




poi 


144 


.458 


.026 


.021 


.037 








Pll 


144 


.001 


.000 


.001 


.001 








P21 


144 


5.046 


1.759 




19. /61** 








P02 


144 


5.506 


195.745*** 


360069.375*** 


378405.077*** 








P12 


144 


.048 


423*** 


.054*** 


.052*** 








P22 


144 


.019 


.256*** 


.053*** 


.052*** 




RE*M 


1 


poi 


72 


5540.631*** 


454326.369*** 


1.038E7*** 


1.557E7*** 








Pll 


72 


.011 


15.816*** 


40.207*** 


40.276*** 








P21 


72 


.014 


15.506*** 


13.080*** 


56.497*** 






2 


P02 


72 


.675 


199.404*** 


716966.573*** 


755012.778*** 








P12 


72 


.201 


.078 


.007** 


.007 








P22 


72 


.116*** 


.529*** 


131*** 


1 34*** 




CR*M 


1 


poi 


72 


.515 


.002 


.002 


.004 








Pll 


72 


.001 


4.520E-5 


6.406E-5 


6.680E-5 








P21 


72 


5.943 


2.889*** 


2.473** 


3.902 






2 


P02 


72 


3.940 


196.384*** 


683436.471*** 


721549.650*** 








P12 


72 


.243 


.407*** 


.004 


.022*** 








P22 


72 


.148*** 


1.340*** 


289*** 


.246*** 




RE*CR*M 


1 


poi 


864 


.348 


.020 


.017 


.030 








Pll 


864 


.001 


.000 


.000 


.000 








P21 


864 


3.917 


1.358 


5.731 


15.895 






2 


P02 


864 


33.150 


884.547*** 


2141981.317*** 


2251208.104*** 








P12 


864 


.072 


.059 


.006 


.005 








P22 


864 


.082 


.433 


.102 


.103 




ERROR 


1 


poi 


1183 


3595.810 


8759.488 


8834975.252 


8871627.167 








(311 


1183 


.245 


11.841 


50.548 


50.564 








P21 


1183 


84.299 


23.465 


28.185 


128.394 






2 


P02 


1 183 


214.134 


135.089 


1384793.426 


1438377.807 








P12 


1183 


16.754 


2.534 


.089 


.150 








P22 


1183 


.659 


1.140 


.197 


.256 




TOTAL 


1 


poi 


2365 


10346.814 


662542.071 


2.575E7 


3.487E7 








(311 


2365 


.288 


64.783 


192.657 


192.913 








P21 


2365 


106.707 


78.742 


85.610 


326.275 






2 


P02 


2365 


261.126 


1804.890 


5756781.519 


6037516.906 








P12 


2365 


18.374 


4.302 


.480 


.365 








P22 


2365 


4.259 


6.818 


1.387 


1.320 



Significant at 0.05 level of significance 



Significant at 0.1 level of significance 
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4.5.1 EFFECT ON p 0 

Consequently in equation 1, it can be inferred that the performances of the estimators are affected by 
autocorrelation under all criteria. The results of the LSD further test visa- vice their estimated marginal means 
revealed that all estimatorsexceptGLS2 are preferred to estimate p 0 at all the levels of autocorrelation. 
In equation 2, the estimators are affected by autocorrelation and correlation between the error terms under all the 
criteria except in bias criterion. The results of the LSD further test visa- vice their estimated marginal means as 
shown in revealed that all estimatorsexceptGLS2 are preferred to estimate Po at all levels of autocorrelation and 
correlation between the error terms. 

EFFECT ON pi 

Consequently in equation 1, it can be inferred that the performances of the estimators are affected by 
autocorrelation under all criteria except for bias. The results of the LSD further test visa- vice their estimated 
marginal means revealed that GLS2 and MLE estimators are preferred to estimate Pi at all the levels of 
autocorrelation. 

In equation 2, the estimators are affected by autocorrelation and correlation between the error terms 
under all the criteria except in bias criterion. The results of the LSD further test visa- vice their estimated 
marginal means revealed that all estimators are preferred to estimate Pi at all levels of autocorrelation and 
correlation between the error terms. 

EFFECT ON p 2 

Consequently in equation 1, it can be inferred that the performances of the estimators are affected by 
autocorrelation and correlation between the error terms. The results of the LSD further test visa- vice their 
estimated marginal means revealed that GLS2 and MLE estimators are preferred to estimate P2 at all the levels 
of autocorrelation. 

In equation 2, the estimators are affected by autocorrelation and correlation between the error terms 
under all the criteria. The results of the LSD further test visa- vice their estimated marginal means revealed that 
SUR and 3SLS estimators are preferred to estimate p 2 at all levels of autocorrelation and correlation between the 
error terms EXCEPT for -0.9 and -0.8 levels of correlation between the error terms under bias that is 
significantly different. 

Summarily, GLS2, MLE, SUR and 3SLS are preferred to estimate the model at sample size of 20 



TABLE4.5.2: ANOVA for the sample size of 30 



n 


sov 


EQN 


Pi 


df 


TYPE III SUM OF SQUARES 












Bias 


Absolute Bias 


Variance 


Mean Square 


30 


RE 


1 


(301 


12 


1368.073*** 


165272.612*** 


1.008E12*** 


1.009E12*** 








(311 


12 


.029*** 


37.228*** 


125.437*** 


127.031*** 








(321 


12 


.075 


29.385*** 


69.897*** 


83.604*** 






2 


(302 


12 


.095 


51.854*** 


6251392.175*** 


6258276.131*** 








(312 


12 


.005 


.334*** 


.025*** 


.025*** 








(322 


12 


.011 


175*** 


.008*** 


.043*** 




CR 


1 


(301 


12 


.147 


10.954 


1.102E8 


1.102E8 








(311 


12 


27 \ 


.043 


.144 


.137 








(321 


12 


1.980*** 


5.721*** 


5.706*** 


5.296*** 






2 


(302 


12 


.200 


13.613*** 


5187309.671*** 


5192414.391*** 








P12 


12 


2.338*** 


.012 


.096*** 


.001 








(322 


12 


.695 


15.467*** 


.018*** 


2.358*** 




M 


1 


(301 


6 


187.891*** 


82932.248*** 


6.302E11 


6.312E11*** 








(311 


6 


.013** 


5.296*** 


7.289*** 


7.403*** 








(321 


6 


.009 


4.316*** 


8.386*** 


8.620*** 






2 


(302 


6 


.007 


9.955*** 


5467535.629*** 


5474905.799*** 








(312 


6 


.029 


.034 


.004*** 


.004 








022 


6 


.001 


.114*** 


.052*** 


.040*** 




RE*CR 




(301 


144 


1.846 


131.775 


1.362E9 


1.363E9 








(311 


144 


.143** 


.163 


1.140 


1.113 








(321 


144 


1.132 


2.945*** 


20.161*** 


18.791*** 








(302 


144 


1.196 


70.924*** 


3.657E7*** 


3.661E7*** 








(312 


144 


.024 


.222 


.016*** 


.016 








(322 


144 


.034 


.108 


.005*** 


.016 
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RE*M 


1 


(301 


72 


7396.149*** 


696530.165*** 


6.041E12*** 


6.050E12*** 








(311 


72 


.012 


17.083*** 


51.754*** 


52.348*** 








(321 


72 


.045 


15.567*** 


46.940*** 


48.988*** 






2 


(302 


72 


583 


43.075*** 


3.735E7*** 


3.739E7*** 








(312 


72 


.002 


.034 


.003** 


.004 








(322 


72 


.028 


.105 


.016*** 


.026 




CR*M 


1 


(301 


72 


.988 


45.595 


6.622E8 


6.623E8 








(311 


72 


203 *** 


.032 


.108 


.103 








(321 


72 


1.491*** 


4.279*** 


4.281 


3.968 






2 


(302 


72 


.863 


40.424*** 


3.189E7*** 


3.193E7*** 








P12 


72 


.012 


.085 


.004*** 


.005 








(322 


72 


.009 


414*** 


.035*** 


.065 




RE*CR*M 


1 


(301 


864 


10.321 


546.890 


8.173E9 


8.174E9 








(311 


864 


.107 


.121 


.855 


.835 








(321 


864 


.853 


2.202 


15.134 


14.102 






2 


(302 


864 


7.135 


249.877*** 


2.190E8*** 


2.192E8*** 








(312 


864 


.001 


.041 


.003 


.003 








(322 


864 


.014 


.174 


.013*** 


.030 




ERROR 


1 


(301 


1183 


3150.131 


4943.259 


4.933E10 


5.003E10 








(311 


1183 


.916 


19.633 


110.203 


112.032 








(321 


1183 


13.579 


13.356 


69.548 


66.923 






2 


(302 


1183 


44.274 


32.901 


2.213E7 


2.212E7 








P12 


1183 


7.545 


6.245 


.041 


.357 








(322 


1183 


123.484 


4.942 


.011 


1.476 




TOTAL 


1 


(301 


2365 


12115.647 


950459.247 


7.739E12 


7.751E12 








(311 


2365 


1.693 


79.598 


296.927 


301.002 








(321 


2365 


19.164 


77.788 


240.059 


250.302 






2 


(302 


2365 


54.357 


512.744 


3.639E8 


3.642E8 








(312 


2365 


9.965 


7.003 


.192 


.413 








(322 


2365 


124.277 


21.562 


.158 


4.059 



4.5.2 EFFECT ON p 0 

Consequently in equation 1, it can be inferred that the performances of the estimators are affected by 
autocorrelation under all criteria. The results of the LSD further test visa- vice their estimated marginal means 
revealed that all estimators are preferred to estimate (3 0 at all the levels of autocorrelation exceptfor GLS2 which 
differed significantly at 0.8, 0.9 and 0.99 autocorrelation levels. 

In equation 2, the estimators are affected by autocorrelation and correlation between the error terms under all the 
criteria except in bias criterion. The results of the LSD further test visa- vice their estimated marginal means 
revealed that all estimators are preferred to estimate (3 0 at all levels of autocorrelation and correlation between 
the error terms exceptfor GLS2 which differed significantly at autocorrelation level of 0.9 and correlation 
between the error terms of 0.99 under bias criterion. 

EFFECT ON pi 

Consequently in equation 1, it can be inferred that the performances of the estimators are affected by 
autocorrelation and correlation between the error terms under all criteria. The results of the LSD further test 
visa- vice their estimated marginal means revealed that GLS2 and MLE estimators are preferred to estimate Pi at 
all the levels of autocorrelation and correlation between the error terms. 

In equation 2, the estimators are affected by autocorrelation and correlation between the error terms under 
variance criterion. The results of the LSD further test visa- vice their estimated marginal means revealed that 
GLS2 and MLE estimators are preferred to estimate Pi at all levels of autocorrelation and correlation between 
the error terms. 

EFFECT ON p 2 

Consequently in equation 1, it can be inferred that the performances of the estimators are affected by 
autocorrelation and correlation between the error terms under all criteria. The results of the LSD further test 
visa- vice their estimated marginal means revealed that GLS2 and MLE estimators are preferred to estimate Pi at 
all the levels of autocorrelation and correlation between the error terms except that we have to be cautious when 
using them at some levels of autocorrelation. 
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In equation 2, the estimators are affected by autocorrelation and correlation between the error terms 
under variance criterion. The results of the LSD further test visa- vice their estimated marginal means revealed 
that all estimators except OLS, GLS2 and MLE estimators are preferred to estimate p 2 at all levels of 
autocorrelation and correlation between the error terms. 



Summarily, GLS2 and MLE estimators are preferred to estimate the model at sample size of 30 

TABLE4.5.3: ANOVA for sample size of 50 



n 


sov 


Pi 


df 


TYPE III SUM OF SQUARES 










Bias 


Abs.Bias 


Var 


MSE 


50 


RE 


(301 


12 


452.571*** 


74575.669*** 


1.764E11*** 


1.770E11*** 






pTl 


12 


.050*** 


18.709*** 


24.791*** 


24.976*** 






P21 


12 


1.014*** 


6.985*** 


1.255*** 


2.662*** 






P02 


12 


.515*** 


35.964*** 


251158.322*** 


252591.912*** 






P12 


12 


4 17*** 


.167*** 


.004*** 


007*** 






P22 


12 


.129** 




.001*** 


.003*** 




CR 


poi 


12 


1.992 


234.178 


1.404E9 


1.406E9 






Pll 


12 


.022*** 




1.780*** 


1.786*** 






P21 


12 


5.131*** 


1 77 


.052 


.030 






P02 


12 


1.353*** 


6.205*** 


161713.711*** 


162539.579*** 






P12 


12 


3.505*** 


1.093*** 


.026*** 


.024*** 






P22 


12 


.221*** 


273*** 


.003*** 


.003*** 




M 


poi 


6 


227.569*** 


24107.884*** 


8.971E10*** 


9.003E10*** 






Pll 


6 


.001 


3.527*** 


2 193*** 


2.209*** 






P21 


6 


.085 


2.311*** 


619*** 


709*** 






P02 


6 


.105*** 


1.285*** 


178487.825*** 


179459.300*** 






P12 


6 


.021*** 


.003 


.002*** 


8.218E-5 






P22 


6 


.023 


307*** 


.010*** 


.012*** 




RE*CR 


poi 


144 


23.036 


2764.733 


1.713E10 


1.714E10 






Pll 


144 


.019 


4.251*** 


12.667*** 


12.699*** 






P21 


144 


1.698 


1.158 


.538*** 


1 049*** 






P02 


144 


2.520*** 


28.684*** 


1365792.064*** 


1373566.001*** 






P12 


144 


.136*** 


.165*** 


.005*** 


.009 






P22 


144 


.058 


.055 


.001 


.001 




RE*M 


poi 


72 


3285.331*** 


280727.544*** 


1.056E12*** 


1.060E12*** 






Pll 


72 


.021 


8.107*** 


10.064*** 


10.139*** 






P21 


72 


.847 


3.449*** 


1.724*** 


1.696*** 






P02 


72 


2.363*** 


10.635*** 


1480512.907*** 


1488990.865*** 






P12 


72 


1 32*** 


2ig*** 


.003*** 


008*** 






P22 


72 


.009 


.064 


.002*** 


.003 




CR*M 


poi 


72 


11.459 


1223.011 


8.561E9 


8.571E9 






Pll 


72 


.016 


.241 


.671 


.673 






P21 


72 


3.879*** 


.469 


.307*** 


.207 






P02 


72 


1.161*** 


6.135*** 


982226.051*** 


987452.884*** 






P12 


72 


.012 


092*** 


.002** 


007*** 






P22 


72 


.010 


.276*** 


.006*** 


007*** 




RE*CR*M 


poi 


864 


137.160 


14657.361 


1.026E11 


1.027E11 






Pll 


864 


.013 


1.423 


4.996 


5.009 






P21 


864 


1.263 


.490 


1.123 


.968 






P02 


864 


13.881*** 


59.010*** 


8150369.132*** 


8196863.794*** 






P12 


864 


.048 


.067 


.002 


.003 






P22 


864 


.008 


.050 


.002 


.002 




ERROR 


poi 


1183 


6296.390 


82375.378 


8.427E11 


8.460E11 






Pll 


1183 


.278 


8.667 


23.218 


23.233 






P21 


1183 


12.234 


8.004 


2.939 


4.219 






P02 


1183 


4.175 


9.073 


1571334.088 


1580748.512 






P12 


1183 


.661 


.922 


.026 


.070 






P22 


1183 


6.136 


1.510 


.012 


.043 




TOTAL 


poi 


2365 


10435.662 


480676.871 


2.295E12 


2.303E12 






Pll 


2365 


.420 


45.888 


80.372 


80.717 






P21 


2365 


26.153 


23.042 


8.557 


11.540 






P02 


2365 


26.079 


157.021 


1.414E7 


1.422E7 






P12 


2365 


4.946 


2.729 


.070 


.128 






P22 


2365 


6.593 


2.811 


.038 


.073 
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4.5.3.1 EFFECT ON p 0 

Consequently in equation 1, it can be inferred that the performances of the estimators are affected by 
autocorrelation under all criteria. The results of the LSD further test visa- vice their estimated marginal means 
revealed that all estimators are preferred to estimate Po at all the levels of autocorrelation exceptfor GLS2 which 
differed significantly at 0.99 autocorrelation levels. 

In equation 2, the estimators are affected by autocorrelation and correlation between the error terms 
under all the criteria. The results of the LSD further test visa- vice their estimated marginal means revealed that 
all estimators are preferred to estimate p 0 at all levels of autocorrelation and correlation between the error terms 
exceptfor GLS2 which differed significantly at autocorrelation level of 0.9 & 0.99 and correlation between the 
error terms of 0.99 under all criteria. 

4.5.3.2 EFFECT ON pi 

Consequently in equation 1, it can be inferred that the performances of the estimators are affected by 
autocorrelation and correlation between the error terms under all criteria. The results of the LSD further test 
visa- vice their estimated marginal means revealed that GLS2 and MLE estimators are preferred to estimate fti at 
all the levels of autocorrelation and correlation between the error terms. 

In equation 2, the estimators are affected by autocorrelation and correlation between the error terms 
under variance criterion. The results of the LSD further test visa- vice their estimated marginal means revealed 
that GLS2 and MLE estimators are preferred to estimate P[ at all levels of autocorrelation and correlation 
between the error terms. 

4.5.3.3 EFFECT ON p 2 

Consequently in equation 1, it can be inferred that the performances of the estimators are affected by 
autocorrelation and correlation between the error terms. The results of the LSD further test visa- vice their 
estimated marginal means revealed that GLS2 and MLE estimators are preferred to estimate p 2 at all the levels 
of autocorrelation and correlation between the error terms. 

In equation 2, the estimators are affected by autocorrelation and correlation between the error terms 
under all the criteria. The results of the LSD further test visa- vice their estimated marginal means revealed that 
all estimators except OLS,GLS2 and MLE estimators are preferred to estimate p 2 at all levels of autocorrelation 
and correlation between the error terms. 



Summarily, GLS2 and MLE estimator are preferred to estimate the model at sample size of 50 

TABLE 4.5.4: ANOVA for sample size of 100 



n 


sov 


Pi 


df 


TYPE III SUM OF SQUARES 










Bias 


Abs.Bias 


Var 


MSE 


100 


RE 


(301 


12 


47.699 


48743.394*** 


7.898E9*** 


7.945E9*** 






(311 


12 


.022*** 


21.108*** 


33.739*** 


33.776*** 






P21 


12 


.014 


7 287*** 


3.435*** 


4.171*** 






(302 


12 


.031 


23.608*** 


27762.719*** 


27802.782*** 






P12 


12 


.004 


.122*** 


.002*** 


.002*** 






(322 


12 


.007 


019*** 


.001*** 


.001*** 




CR 


(301 


12 


.001 


.001 


.000 


.000 






(311 


12 


QH *** 


.002 


.006 


.006 






P21 


12 


1.366*** 


1.539*** 


.352*** 


320*** 






P02 


12 


.486 


139*** 


12228.857*** 


12185.714*** 






(312 


12 




.692*** 


.018*** 


.016*** 






|322 


12 


.057 


1.315*** 


.002*** 


.019*** 




M 


(301 


6 


12.616 


13036.510*** 


3.909E9*** 


3.932E9*** 






(311 


6 


.005*** 


3.739*** 


2.145*** 


2.151*** 






P21 


6 


.000 


1.601*** 


.568*** 


.585*** 






(302 


6 


.044 


.218*** 


13231.759*** 


13252.284*** 






P12 


6 


.003 


.058*** 


.002*** 


.002*** 






P22 


6 


.002 


.095*** 


.004*** 


.005*** 




RE*CR 


poi 


144 


.001 


.008 


.002 


.002 






Pll 


144 


.022*** 


.026 


.077 


.077 






P21 


144 


.810*** 


.784*** 


1.083*** 


97g*** 






P02 


144 


.204 


15.466*** 


147798.362*** 


148004.003*** 






P12 


144 


.011 


.053*** 


.002*** 


.002*** 






P22 


144 


.002 


.021 


.000*** 


.000 




RE*M 


poi 


72 


151.558 


156345.485*** 


4.691E10*** 


4.719E10*** 






Pll 


72 


.015*** 


9.522*** 


13.766*** 


13.781*** 
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(321 


72 


.005 


4.081*** 


2.357*** 


2 457*** 






(302 


72 


.167 


2.498*** 


158513.192*** 


158738.340*** 






P12 


72 


.011 


.044*** 


.001*** 


.001*** 






(322 


72 


.009 


029*** 


.001*** 


.001*** 




CR*M 


(301 


72 


.001 


7.294E-6 


5.759E-5 


6.569E-5 






(311 


72 


.009*** 


.001 


.004 


.004 






(321 


72 


1.024*** 


1.158*** 


.264*** 


.240*** 






|302 


72 


.063 


1 170*** 


73197.801*** 


73298.063*** 






(312 


72 


.001 


.052*** 


.001*** 


.001*** 






(322 


72 


.002 


123*** 


.002*** 


.002*** 




RE*CR*M 


(301 


864 


.001 


.000 


.001 


.001 






(311 


864 


.018 


.016 


.058 


.057 






(321 


864 


.610 


.587 


.813 


.734 






(302 


864 


.899 


13.854*** 


877548.891*** 


878767.713*** 






P12 


864 


.011 


.029 


.001 


.001 






P22 


864 


.005 


.038 


.001*** 


.001 




ERROR 


poi 


1183 


3458.358 


1982.236 


5.755E10 


5.739E10 






Pll 


1183 


.034 


8.157 


27.089 


27.121 






(321 


1183 


1.548 


1.645 


1.959 


2.328 






P02 


1183 


28.879 


4.729 


142202.156 


142616.134 






P12 


1183 


.422 


.148 


.005 


.005 






P22 


1183 


13.178 


.257 


.000 


.009 




TOTAL 


poi 


2365 


3670.242 


220115.863 


1.163E11 


1.165E11 






(311 


2365 


.137 


42.568 


76.883 


76.972 






P21 


2365 


5.378 


18.685 


10.833 


11.814 






P02 


2365 


30.774 


61.678 


1452558.902 


1454739.728 






P12 


2365 


.510 


1.199 


.031 


.030 






P22 


2365 


13.260 


1.905 


.012 


.039 



4.5.4.1 EFFECT ON p 0 

Consequently in equation 1, it can be inferred that the performances of the estimators areaffected by 
autocorrelation under absolute bias, variance and mean square error criteria. The results of the LSD further test 
visa- vice their estimated marginal means revealed that all estimators are preferred to estimate (3o at all the levels 
of autocorrelation exceptfor GLS2 which differed significantly at 0.99 autocorrelation levels. 
In equation 2, the estimators are affected by autocorrelation and correlation between the error terms under all the 
criteria. The results of the LSD further test visa- vice their estimated marginal means revealed that all estimators 
are preferred to estimate p 0 at all levels of autocorrelation and correlation between the error terms exceptfor 
GLS2 which differed significantly at autocorrelation level of 0.99 and correlation between the error terms of - 
0.99 and +0.99 under all criteria considered. 

4.5.4.2 EFFECT ON p, 

Consequently in equation 1, it can be inferred that the performances of the estimators areaffected by 
autocorrelation and correlation between the error terms under all criteria. The results of the LSD further test 
visa- vice their estimated marginal means revealed that GLS2 and MLE estimators are preferred to estimate Pi at 
all the levels of autocorrelation and correlation between the error terms. 

In equation 2, the estimators are affected by autocorrelation and correlation between the error terms 
under all the criteria. The results of the LSD further test visa- vice their estimated marginal means revealed that 
GLS2 and MLE estimators are preferred to estimate P[ at all levels of autocorrelation and correlation between 
the error terms. 

4.5.4.3 EFFECT ON p 2 

Consequently in equation 1, it can be inferred that the performances of the estimators are affected by 
autocorrelation and correlation between the error terms under all criteria. The results of the LSD further test 
visa- vice their estimated marginal means revealed that GLS2 and MLE estimators are preferred to estimate P2 at 
all the levels of autocorrelation and correlation between the error terms, although they too are significantly 
different at some limited levels of autocorrelation. 

In equation 2, the estimators are affected by autocorrelation and correlation between the error terms 
under variance criterion. The results of the LSD further test visa- vice their estimated marginal means as shown 
in Table EB4.8 of appendix 5 revealed that all estimators except OLS, GLS2 and MLE estimators are preferred 
to estimate P2 at all levels of autocorrelation and correlation between the error terms. 
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Summarily, GLS2, SUR and MLE estimator are preferred to estimate the model at sample size of 100 
TABLE4.5.5: ANOVA for sample size of 250 



n 


sov 


Pi 


df 


TYPE III SUM OF SQUARES 










Bias 


Abs.Bias 


Var 


MSE 


250 


RE 


(301 


12 


1.478*** 


4632.931*** 


1.059E8*** 


1.059E8*** 






(311 


12 


.030*** 


6.158*** 


2.703*** 


2.730*** 






P21 


12 


.001 


2.812*** 


.540*** 


.658*** 






(302 


12 


.003 


6.709*** 


83.297*** 


83.466*** 






P12 


12 


.319 


.035*** 


17.412 


.001** 






(322 


12 


.002 


021 *** 


8.986E-5*** 


.000*** 




CR 


(301 


12 


.008 


6.104E-5 


8.877E-5 


8.761E-5 






(311 


12 


.001 


5.698E-5 


3.105E-5 


3.798E-5 






(321 


12 


.205*** 


.627*** 


.062*** 


.058*** 






(302 


12 


.338*** 


295*** 


30.040 


29.626 






(312 


12 


.778** 


.356*** 


17.229 


.006*** 






(322 


12 


.036 


303*** 


.000*** 


.001*** 




M 


(301 


6 


.190*** 


873.346*** 


5.133E7*** 


5.133E7*** 






(311 


6 


7.732E-5 


1.168*** 


229*** 


230*** 






P21 


6 


.001 


.583*** 


102*** 


.105*** 






P02 


6 


.001 


.021 


28.674 


28.775 






(312 


6 


.182 


.007 


8.763 


.000 






(322 


6 


.001 


073*** 


.001*** 


.001*** 




RE*CR 


(301 


144 


.005 


.000 


.001 


.001 






(311 


144 


.012 


9.808E-5 


.000 


.000 






(321 


144 


.176 


.311** 


.155*** 


.148*** 






P02 


144 


.027 


2.815*** 


352.222 


353.150 






P12 


144 


3.609 


.015 


209.071 


.000 






P22 


144 


.005 


.012 


6.516E-5*** 


.000 




RE*M 


poi 


72 


2.434 


10508.921*** 


6.160E8*** 


6.161E8*** 






pil 


72 


on*** 


2.858*** 


1.109*** 


1.120*** 






P21 


72 


.001*** 


1.528*** 


.355*** 


.376*** 






P02 


72 


.003 


.195 


344.710*** 


345.580*** 






P12 


72 


1.811 


.071 


104.517 


.001 






P22 


72 


.002 


042*** 


.000*** 


.000*** 




CR*M 


poi 


72 


.006 


4.58 1E-5 


6.662E-5 


6.574E-5 






Pll 


72 


.001 


3.947E-5 


2.313E-5 


2.817E-5 






P21 


72 


.152** 


.469*** 


.046 


.044 






P02 


72 


.005 


.090 


162.122 


162.533 






P12 


72 


1.763 


.013 


104.616 


.000 






P22 


72 


.006 


075*** 


.001*** 


.001*** 




RE*CR*M 


poi 


864 


.004 


.000 


.001 


.001 






Pll 


864 


.009 


8.591E-5 


.000 


.000 






P21 


864 


.131 


.233 


.116 


.111 






P02 


864 


.055 


1.121 


1945.383 


1950.633 






P12 


864 


21.396 


.030 


1254.431 


.001 






P22 


864 


.004 


.021 


.000*** 


.000 




ERROR 


poi 


1183 


3.854 


16061.289 


7.821E8 


7.822E8 






Pll 


1183 


.082 


3.122 


2.240 


2.268 






P21 


1183 


1.781 


2.054 


.749 


.792 






P02 


1183 


3.320 


7.305 


2945.654 


2953.427 






P12 


1183 


37.560 


1.475 


1717.611 


.034 






P22 


1183 


3.873 


.314 


6.234E-5 


.003 




TOTAL 


poi 


2365 


7.981 


32077.122 


1.555E9 


1.556E9 






Pll 


2365 


.146 


13.305 


6.281 


6.348 






P21 


2365 


2.449 


8.618 


2.125 


2.292 






P02 


2365 


3.756 


18.552 


5892.407 


5907.490 






P12 


2365 


67.419 


2.003 


3433.660 


.042 






P22 


2365 


3.929 


.863 


.003 


.007 
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4.5.5.1 EFFECT ON p 0 

Consequently in equation 1, it can be inferred that the performances of the estimators areaffected by 
autocorrelation under all criteria. The results of the LSD further test visa- vice their estimated marginal means 
revealed that all estimators are preferred to estimate Po at all the levels of autocorrelation exceptfor GLS2 which 
differed significantly at 0.99 autocorrelation levels. 

In equation 2, the estimators are affected by autocorrelation under variance and mean square error 
criteria. The results of the LSD further test visa- vice their estimated marginal means revealed that all estimators 
are preferred to estimate p 0 at all levels of autocorrelation exceptfor GLS2 which differed significantly at 
autocorrelation level of 0.99 in both criteria considered. 

Summarily, we can infer that all the estimators are preferred to estimate p 0 except GLS2 at all five sample sizes 
under consideration. 

4.5.5.2 EFFECT ON pi 

Consequently in equation 1, it can be inferred that the performances of the estimators areaffected by 
autocorrelation under all criteria. The results of the LSD further test visa- vice their estimated marginal means 
revealed that GLS2 and MLE estimators are preferred to estimate P[ at all the levels of autocorrelation. 

In equation 2, the estimators are neither affected by autocorrelation nor correlation between the error 
terms under all criteria. 

Summarily, we can infer that GLS2 and MLE estimators are preferred to estimate P [ at all five sample 
sizes under consideration and at all levels of autocorrelation and correlation between the error terms. 

4.5.5.3 EFFECT ON p 2 

Consequently in equation 1, it can be inferred that the performances of the estimators areaffected by 
autocorrelation and correlation between the error terms under all criteria. The results of the LSD further test 
visa- vice their estimated marginal means revealed that GLS2 and MLE estimators are preferred to estimate p 2 at 
all the levels except at -0.99 and +0.99 correlation between the error terms under absolute bias. 

In equation 2, the estimators are affected by autocorrelation and correlation between the error terms 
under variance criterion. The results of the LSD further test visa- vice their estimated marginal means as shown 
in revealed that all estimators except OLS, GLS2 and MLE estimators are preferred to estimate p 2 at all levels of 
autocorrelation and correlation between the error terms. 

Summarily, we can infer that GLS2 and MLE estimators are preferred to estimate p 2 . 

Summarily, MLE estimator is preferred to estimate the model at sample size of 250 

Conclusively, the estimator of MLE is preferred to estimate all the parameters of the model in the presence of 
correlation between the error terms and autocorrelation at all the sample sizes. 

Figure 7: Performaces of the estimators using M5E(B L1 } at different levels 
of sanple size, correlation bet the error term and autocorrelation at 

CR=-0.99 




— ♦ — OLS 
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— * — MLE 

-M-MLILTIREG 

-*HFIML 

-•-5LJR 

—I — 3SLS 
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In fig. 7, the plot of MSE values against different sample sizes for all the estimators revealed appreciable 
increase in efficiency (lower MSE) of the estimators as sample size increases with MLE estimator showing a 
better performance over GLS2. 

V. SUMMARY OF THE FINDINGS 
5.1.5 When there is correlation between the error terms and Autocorrelation 

The summary of results from the Analysis of variance tables of the criteria showing the performances of the 
estimators and sample sizes on parameters of the two equation model when there is presence of correlation 
between the error terms and autocorrelation are presented in Table 5.1.5 

Table 5.1.5: Summary of results when there is presence of correlationbetweenthe error terms and 



autocorrelation 



n 


EQ 

N 


PARAM 
ETERS 


PREFERRED 


OVERALL ASSESSMENT 


MOST 

PREFERRED 


20 


1 


Poi 


All except CORC 


CORC, MLE 


MLE 


Pn 


CORC, MLE 


Pai 


CORC, MLE 


2 


P 02 


All except CORC 


All except CORC 


P 12 


All 


P22 


SUR, 3SLS 


30 


1 


Poi 


All except CORC 


CORC, MLE 


MLE 


Pll 


CORC, MLE 


P21 


CORC, MLE 


2 


P02 


All except CORC 


MLE.SUR,3SLS 


P 12 


CORC, MLE 


P22 


MulReg,FIML,SUR,3SLS 


50 


1 


P01 


All 


CORC, MLE 


CORC, MLE 


Pll 


CORC, MLE 


P21 


CORC, MLE 


2 


P02 


All 


CORC,MLE.SUR,3SLS 


Pl2 


CORC, MLE 


P22 


MulReg,FIML,SUR,3SLS 


100 


1 


P01 


All except CORC 


CORC, MLE 


MLE 


Pll 


CORC, MLE 


P21 


CORC, MLE 


2 


P02 


All except CORC 


MLE,MulReg,FIML,SUR,3SLS 


P 12 


CORC, MLE 


P22 


MulReg,FIML,SUR,3SLS 


250 


1 


P01 


All except CORC 


CORC, MLE 


MLE 


Pll 


CORC, MLE 


P21 


CORC, MLE 


2 


P02 


All except CORC 


All except CORC 


P 12 


All 


P22 


MulReg,FIML,SUR,3SLS 



From table 5.1.5 when there is presence of correlation between the error terms and autocorrelation in 
the model under the equation 1 in all the five sample sizes, all the estimating methods except CORC are equally 
good in estimating the parameters P 0 i, meanwhile, for parameters P nandp 21 CORC and MLE estimators are 
good for their estimation thus it can be concluded that MLE estimating method is preferred in estimating all the 
model parameters in equation 1 . 

Under equation 2, it was observed that all estimation methods except CORC are good in estimating all 
the parameters of the model at all level of the sample sizes. 
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However, observing the two equations together, we can conclude that MLE is the most preferred in 
estimating all the parameters of the two equations among all the estimation methods used. 

5.2. RECOMMENDATION 

The research work has revealed that MLE method of estimation is the most preferred estimator in 
estimating all the parameters of the model based on the four criteria used namely; Bias, Absolute Bias, Variance 
and Mean Square Error under the five level of sample sizes considered. It can therefore be recommended that 
when the validity of other correlation assumptions cannot be authenticated in seemingly unrelated regression 
model, the most preferred estimator to use is MLE. Meanwhile, for any SUR model without any form of 
correlation, SUR estimation method is most preferred. 

5.3. SUGESTION FOR FURTHER STUDY 

This study considered two- equation model with two depended variables in each equation, a future research may 
consider situation in which more than two equations and as many depended variables as possible. 

One may still consider a Bayesian estimation approach as one of the estimation methods in order to test its own 
potential. 
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